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Abstract 

Paul Van Dooren [19j constructed an algorithm for the compu- 
tation of all irregular summands in Kronecker's canonical form of a 
matrix pencil. The algorithm is numerically stable since it uses only 
unitary transformations. 

We construct a unitary algorithm for computation of the canonical 
form of the matrices of a chain of linear mappings 

Vi V2 Vt 

and extend Paul Van Dooren's algorithm to the matrices of a cycle of 
linear mappings 




where all Vi are complex vector spaces and each line denotes — > or 
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1 Introduction 



All matrices and vector spaces are considered over the field C of complex 
numbers. 

By the theorem on pencils of matrices (see Sect. XII]), every pair of 
p X q matrices reduces by transformations of simultaneous equivalence 



{A,, A2) ^ {R-^AiS, R-^A^S) 



(1) 



{R and S are arbitrary nonsingular matrices) to a direct sum, determined 
uniquely up to permutation of summands, of pairs of the form 



(/„,J„(A)), (J„(0),/„), (F„,G„), (F^,^^), 



(2) 



where 



1 






1 



1 






1 



n > 1, 



(3) 
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0" 




'B 


0" 


( 





C 







D 



are (n — 1) x n matrices, and Jn(A) is a Jordan block. The direct sum of 
pairs is defined by 



(A, B) © (C, D) = {A®C,B®D) 



Note that Fi and Gi in ([3]) have size 0x1. It is agreed that there exists 
exactly one matrix, denoted by 0„o; of size n x and there exists exactly one 
matrix, denoted by Oon, of size x n for every nonnegative integer n; they 
represent the linear mappings ^ C" and C" ^ and are considered as 
zero matrices. Then 



pq 







mO 



'M„ " 






Opo" 






_ 0^0. 






OmO 




Orraq 



and 



Mpq © Oo„ 



Mp, 



0, 



On 



Oog Oon 



[Mp, Op„] 



for every p x q matrix M 



pq- 
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p. Van Dooren [12] constructed an algorithm that for every pair (A, B) 
of p X q matrices calculates a simultaneously equivalent pair 



{Ai,Bi)®---®{Ar,Br)®{C, D), 

where all {Ai, Bj) are of the form 

(4,J„(0)), (J„,(0),/„), iF^,Gl), 

and the matrices C and D are nonsingular. The pair (C, D) is called a regular 
part of {A, B) and is simultaneously equivalent to a direct sum of pairs of 
the form J„(A)) with A 7^ 0. This algorithm uses only transformations 

with unitary R and S, which is important for its numerical stability. 

In this article we construct a unitary algorithm for computation of the 
canonical form of the system of matrices of a chain of linear mappings 

Vi V2 • • ■ Vt (4) 

(see Proposition 14.11) and extend Van Dooren's algorithm to the matrices of 
a cycle of linear mappings 

At-2 At-l 

A: V . V2 Ft-i ^ Vt , t ^ 2, 




(see Theorem 16. ip . where each line is the arrow — > or the arrow < — and 
Vi, . . . ,Vt are vector spaces. 

For instance, the linear mappings Ai and A2 of a cycle 

A, 

Vi-. :V2 
A2 

are represented by a pair of matrices {Ai,A2) with respect to bases in Vi 
and V2, and a change of the bases reduces this pair by transformations of 
simultaneous equivalence ([1]); in this case our algorithm coincides with Van 
Dooren's algorithm. 

Similarly, the linear mappings Ai and A2 of a cycle 

Ai 

Vi , ' V2 

A2 
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are represented by a pair (^1,^42), and a change of the bases in Vi and V2 
reduces this pair by transformations of contragredient equivalence 



The direct sum of the cycle ([5]) and a cycle 




with the same orientation of arrows is the cycle A ® A': 



Ai © A[ 



V2 © V' 



Ao © ^9 



At-i © Al, 



t-l 



At®A!^ 



A cycle A of the form ([5]) is called regular if all Ai are bijections; other- 
wise it is called singular. By a regularizing decomposition of A, we mean a 
decomposition 

A = v®---®g®v, (6) 

where X>, . . . , ^ are direct-sum-indecomposable singular cycles and P is a 
regular cycle. 

In Section [2] we recall notions of quiver representations; they allow to 
formulate our algorithms pictorially. 

In Section [3] we recall the classification of chains (jlj) and cycles ([5]) of 
linear mappings. The classification of cycles of linear mappings was obtained 
by Nazarova [15j and, independently, by Donovan and Freislich [5j (see also 
[7], Theorem 11.1). 

In Section m we construct an algorithm that gets the canonical form of the 
matrices of a chain of linear mappings using only unitary transformations. 

In Sections and we construct an algorithm that gets a regularizing 
decomposition (jS]) of a cycle of linear mappings using only unitary transfor- 
mation^. The singular summands P, . . . , ^ will be obtained in canonical 
form. 

^This improves the numerical stability of the algorithms. Nevertheless, this does not 
guarantee that the computed structure of the cycle coincides with its original structure. 
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The canonical form of the (nonsingular) matrices Pi, Pt of the regular 
summand 



V: 




Vt 

in ([6]) is not determined by this algorithm. We may compute it as follows. 
We first reduce Pi to the identity matrix changing the basis in the space U2- 
Then we reduce P2 to the identity matrix changing the basis in the space t/3, 
and so on until obtain 

Pi = ... = P^_i = J„. (7) 

At last, changing the bases of all spaces t/i, . . . , [/^ by the same transition ma- 
trix S (this preserves the matrices ([7])), we can reduce the remaining matrix 
Pt to a nonsingular Jordan canonical matrix $ by similarity transformations 
S^^PtS. Clearly, the obtained sequence 

(4, ...,/„,$) 

is the canonical form of the matrices of V. 

2 Terminology of quiver representations 

The notion of a quiver and its representations was introduced by Gabriel [6] 
(see also [3, Section 7]) and admits to formulate classification problems for 
systems of linear mappings. A quiver is a directed graph; loops and multiple 
arrows are allowed. Its representation A over C is given by assigning to each 
vertex v a complex vector space and to each arrow a : u ^ v a. linear 
mapping Aa : — ^ K of the corresponding vector spaces. 
For instance, a representation of the quiver 



is a system of linear mappings 




The number 

dim^, A := dim Vy 

is called the dimension of A at the vertex v, the set of these numbers 

dim^ := {dimK;}^ 

is called the dimension of A. 

Two representations A and A' are called isomorphic if there exists a set 
S of linear bijections Sy : Ay ^ A^ (assigned to all vertices v) transforming 
A to A'. That is, the diagram 



Vy 



Sy 



(8) 



A' 

V' — ^ v 

U V 



must be commutative {A!^Su = SyAa) for every arrow a : u — ^ v. In this 
case we write 

S = {Sy} :A^A' and A- A'. (9) 

The direct sum of A and A' is the representation A® A' formed by Vy © V^ 
and Aa®A'^. 

The following theorem is a well-known corollary of the Krull-Schmidt 
theorem P, Theorem 1.3.6] and holds for representations over an arbitrary 
field. 

Theorem 2.1. Every representation of a quiver decomposes into a direct sum 
of indecomposable representations uniquely, up to isomorphism of summands. 

Every representation of a quiver over C is isomorphic to a representation, 
in which the vector spaces Vy assigned to the vertices all have the form 
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C © ■ ■ ■ © C. Such a representation of dimension {dy} with dy E {0,1,2, . . .} 
is called a matrix representatioi^ and is given by a set A of matrices G 
(^dyxdu assigned to the arrows a : u — > v. We will consider mainly matrix 
representations. 

For every matrix representation A = {A^} of a quiver Q, we define the 
transpose matrix representation 

A^ = {Al} (10) 

of the quiver obtained from Q by changing the direction of each arrow. 
Clearly, 

§ = {Sy} : A ^ 1 implies §^ = {S^} : ^ A^ (11) 

The systems of linear mappings (jlj) and may be considered as repre- 
sentations of the quivers 

£: 1^2^ ... ^(t-l)^t (12) 




with the same orientations of arrows as in (jl]) and ([5]). The quiver (fT3l) will 
be called a cycle; the symbol C will always denote the cycle (IT^ . 

If A is a matrix representation of a quiver with an indexed set of arrows 
{a, I i G /}, we will write Ai instead of Aq,-. So a matrix representation A of 
the cycle C is given by a sequence of matrices 

A= (Ai,...,A). 

matrix representation also arises when we fix bases in all the spaces of a represen- 
tation. As follows from ([8]), two matrix representations are isomorphic if and only if they 
give the same representation but in possible different bases. 



7 



3 Classification theorems 



In this section, we recall the classification of representations of the quivers 
(fT2|) and (fT3l) . and mention articles considering special cases. Some of these 
articles are little known outside of representation theory. 

We first consider the cycles of length 2. The representations of the cycle 
1 =^ 2 were classified by Kronecker [12j in 1890 (see also [HJ Sect. V] or [TJ 
Sect. 1.8]): every pair of p x g matrices is simultaneously equivalent to a 
direct sum of pairs of the form ([2]) . A simple and short proof of this result 
was obtained by Nazarova and Roiter [T6] . 

A classification of representations of the cycle 1^2 was obtained by 
Dobrovol'skaya and Ponomarev [3] in 1965: every matrix representation is 
isomorphic to a direct sum, determined uniquely up to permutation of sum- 
mands, of matrix representations of the form 

(4,J„(A)), (J„(0),4), (Fr^Gl), {F^,Gn) (14) 

(see ([3])). Over an arbitrary field, the Jordan block Jn(A) is replaced by a 
Frobenius block 

1 



where 

x" + aix^''^ + ■ ■ ■ + a„_ix + a„ = p(x)* 

for some irreducible polynomial p{x) and some integer t. This result was 
proved again by Rubio and Gelonch [T7] in 1992, Olga Holtz [TU] in 2000, and 
Horn and Merino [11] in 1995; the last article also contains many applications 
of this classification. 

A classification of systems of linear mappings of the form 

i T 

was given by Nazarova [13J in 1961 over the field with two elements, and by 
Nazarova [H] in 1967 over an arbitrary field. 
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A quiver is said to be of tame type if the problem of classifying its rep- 
resentations does not contain the problem of classifying pairs of matrices up 
to simultaneous similarity. If a quiver Q is not of tame type, then a full 
classification of its representations is impossible since it must contain a clas- 
sification of representations of all quivers, see [HI Sect. 3.1] or [3], Sect. 2]. 
Nevertheless, each particular representation of Q can be reduced to canonical 
form, see [2] or [13, Sect. 1.4]. 

Nazarova [T5] and, independently, Donovan and Freislich [5] in 1973 clas- 
sified representations of all quivers of tame type (see also [H Sect. 11]). In 
particular, they classified representations of the cycle (fT3|l . which is of tame 
type (see this classification also in [TJ Theorem 11.1]). This classification 
is not mentioned in many articles on linear algebra and system theory that 
study its special cases (for instance, in the article by Gelonch [9] contain- 
ing the classification of representations of the cycle f|T3l) with orientation 
1^2^ ,t^l). 

Gabriel [6] (see also [3, Sect. 11]) classified representations of all quivers 
having a finite number of nonisomorphic indecomposable representations. In 
particular, he classified representations of the quiver f[T^ . 

Now we formulate theorems that classify representations of the quivers 
and (US]). 

For every pair of integers such that 1 ^ i ^ j ^ t, we define the 
matrix representation 

I 1 . Ji /i . 

Lij : 1 • • • I ■ ■ ■ j ■ ■ ■ t (15) 

of dimension (0, . . . , 0, 1, . . . , 1, 0, . . . 0) of the quiver (fT2!) . By the next theo- 
rem, which holds over an arbitrary field, the representations hij form a full 
set of nonisomorphic indecomposable matrix representations of (JT2l) . 



Theorem 3.1 (see [6J). For every system of linear mappings there 
are bases of the spaces Vi, . . . ,Vt, in which the sequence of matrices of 
Ai, . . . , At-i is a direct sum of sequences (0, . . . , 0, Ji, . . . , Ji, 0, . . . , 0) of di- 
mension (0, . . . , 0, 1, . . . , 1, 0, . . . 0). This sum is determined by the system 
(jl]) uniquely up to permutation of summands. 

The classification of representations of a cycle f|T3l) follows from Theorem 
12.11 and the next fact: if a matrix representation of this cycle is direct-sum- 
indecomposable, then at least t — 2 of its matrices are nonsingular. Clearly, 
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these t — 2 matrices reduce to the identity matrices and the remaining two 
matrices reduce to the form ([2]) or fll4p depending on the orientation of their 
arrows. This gives the following theorem. 

Theorem 3.2 (see [S] or [U])- For every system of linear mappings 
there are bases in the spaces Vi, . . . ,Vt, in which the sequence of matrices of 
Ai, . . . , At is a direct sum, determined by ([5]) uniquely up to permutation of 
summands, of sequences of the following form {the points denote sequences 
of identity matrices or Oqo) 

(i) (J„(A),...) withX^Q; 

(ii) (. . . , Jri(O), . . .) with J„,(0) at the place i G {1, . . . , t}; 

(iii) (. . . , v4j, . . . , v4j, . . where Ai and Aj depend on the direction of the 
mappings Ai and Aj in the sequence 

{see dS])) as follows: 

^ \{Fn,Gn) or{F^,G'^) if Ai and Aj have opposite directions, 
^ " ''^\{F^,Gl) or{F^,Gr,) otherwise. 

This theorem, with a nonsingular Frobenius block instead of Jn(A) in (i), 
holds over an arbitrary field. 

In the remaining part of this section, we recall Gabriel and Roiter's con- 
struction [3, Sect. 11.1] of summands (ii) and (iii). 

For every integer n, denote by [n] the natural number such that 

1 ^ [n] ^ t and [n] = n mod t . 

Let 

/ — (/ + 1) — {1 + 2) — r, 1^/^t, (16) 

be a "clockwise walk" on the cycle f|T3|) ) that starts at the vertex I, passes 
through the vertices 

[/ + 1], [/ + 2], ... , [r-1], 
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and stops at the vertex [r]. This walk determines the representation A of C 
in which each space K is spanned hj alii E {1,1 + 1, ... ,r} such that [i] = v: 

Vv = ^ i ^ r, [i] = v), 

and all the nonzero actions of linear mappings Aai , ■ ■ ■ , Aat on the basis 
vectors are given by ffTBl) . The matrices of Aa^, • • • , Aat in these bases form 
a matrix representation denoted by 



G, 



Ir- 



(17) 



Example 3.1. The walk 



1 
7 



3, 
9 



on the cycle 



C: 1 



2 
6 



3 
5 



determines the representation 



^1,9: (1,7) 

[0^ 



(2,8) 



(6) 



(3,9) 



(5) 



(4) 



Lemma 3.1 (see [TJ Sect. 11.1]). The set of all Gi^ coincides with the set of 
matrix representations of the form (ii) and (iii).' 

(a) G;^. with r ^ I — 1 mod t is the matrix representation (iii) of dimension 
{di, . . . , dt) , where di is the number of n G {1,1 + 1, . . . ,r} such that 
[n] = i. {Note that all representations of the form (iii) have distinct 
dimensions and so they are determined by their dimensions.) 
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(b) Gi^i^i+pt = {Ip,..., Ip, Jp{0),lp, ...,Ip), where Jp(0) is at the [I - l]-st 
place. 



We will use the following notation. If all arrows in a representation 



U2 



v4i 



have the same orientation, then instead of (1181) we will write 



where 



A= I 



Ui 
Un 



[Ai\. . .\An] if Ml > V, U2 > V,...,Un > V, 

if Ui < V, U2 < V, . . . ,Un * V. 



An 



(18) 



(19) 



(20) 



The partition of A into strips is fully determined by the dimensions of (|T8 
at the vertices ui, . . . , m„. 



4 Chains of linear mappings 

In this section we give an algorithm that calculates the canonical form of the 
matrices of a chain of linear mappings (jlj) using only unitary transformations. 

Choosing bases in the spaces Vi, . . . , we may represent a system of 
linear mappings (jl]) by the sequence of matrices A = (Ai, . . . , At-i). We will 
consider this sequence as a matrix representation 

At-i , 



Ai A2 
1 — - 2 — - 



21 



of the quiver f|T2l) . 

For every vertex i, a change of the basis in Vi changes A. This transfor- 
mation of A will be called a transformation at the vertex i. It will be called a 
unitary transformation if the transition matrix to a new basis of Vi is unitary. 
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The algorithm for chains: 

Let A be a matrix representation fl2T|) of dimension 

dim A = {di, . . . ,dt) 

of the quiver f|T2|) . We will sequentially split A into representations of the 
form f[T^ . 

Step 1: By unitary transformations at vertices 1 and 2, we reduce Ai to the 
form 

To H] 







(22) 



where if is a nonsingular k x k matrix. These transformations change A2] 
denote the new matrix by Ag. Denote also by Vi the set consisting of di — k 
representations of the form Ln. 

Next we will transform the representation A into a representation Mi of 
a "split" quiver depending on the direction of ai in f[T^ as follows: 

Case ai : 1 — > 2. Then 



{k copies) 



1^2 
'l'^'2' 



A' 



3 ^3 4 ^4 



((i2 — copies) 

(there are k fragments of the form 1 — >2 3 and d2 — k fragments of the 

form 2 — 3). The direction of the arrows is the same as in the quiver ( fT2l) . 

Case ai : 1 < — 2. Then 



(d2 — k copies) 



A' 



{k copies) 



1^2 
1^2 



3 ^3 4 ^4 



Step r (1 < r < t): Assume we have constructed in the step r — 1 the set 
Vr-i consisting of representations of the form Lj^, 1 ^ i ^ j < r, and a 
quiver representation M,,_i: 
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(ki copies) pi (pi + 1) ■ — — r n 
{k2 copies) p2 — {p2 + 1) — — r ~~~~ 

{kr copies) pr {pr + 1) — ^ " " " — ^ r 
in which every 

h h 
Pi {Pi + 1) 



(r + 11 



t (23) 



repeats k^ times, fci + ■ ■ ■ + A;^ = dr, all A;^ ^ 0, and 

{pi,P2,---,Pr] = {1,2, ...,r}. 
The direction of the arrows is the same as in the quiver (|T2 



Case ar : r — > r + 1 (see (fT2l)). We divide into r vertical strips of 
sizes ki, k2, ■ ■ ■ , kr and reduce ^4' to the form 



(24) 








Hi 


* 


* 


* 


* 




* 


* 













H2 


* 


* 




* 


* 



















H3 




* 


* 


























































(where each Hi is a nonsingular /j x 1^ matrix and each * is an unspecified ma- 
trix) starting from the first vertical strip by unitary column-transformations 
within vertical strips and by unitary row-transformations. These row- 
transformations are transformations at the vertex r + 1 and they change 
Ar+i] denote the obtained matrix by A'^^^. Denote also by Vr the set ob- 
tained from Vr-i by including k^ — li representations of the form Lp.^ for all 
i = 1, . . . , r. Construct the quiver representation 



{li copies) pi 



h 



h 



(r + 1) 



V+l 



h 



2) 



lr + 2 



{Ir copies) Pr 

{dr+1 — h — ■ ■ ■ — Ir copies) 



(r + 1) 
(r + 1) 



(Hence, k^ — li representations 



•piT 



Pi 
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for each i = 1, . . . , r "break away" from the representation (1231) and join to 
the set Vr-i-) In particular, if r = t — 1, then takes the form 



(/i copies) pi 



h 



h 



t 



h 



{k-i copies) pt-i— —t 

{dt — h — ■ ■ ■ — k-i copies) t 



(25) 



Case Ur : r < — r + 1. We partition A'^ into r horizontal strips of sizes 
ki, k2, . . . , K and reduce A' to the form 



Br 



Hi * 
0* 









■■ 


Hr-2 


* 


* 








■■ 





* 


* 








•■ 





Hr-1 


* 








■■ 








* 








•• 








Hr 








•• 












(26) 



(where each Hi is a nonsingular li x li matrix) starting from the lower strip, 
by unitary row-transformations within horizontal strips and by unitary 
column-transformations. These column-transformations are transformations 
at the vertex r + 1 and they change Ar+i, denote the obtained matrix by 
Ar_^_i. Denote also by Vr the set consisting of the elements of Vr-i and 
ki — li representations of the form Lp.^ for all i = 1, . . . ,r. Construct the 
quiver representation 

{dr+i — li — ■ ■ ■ — Ir copies) (r + 1) 
(/i copies) pi — — ■ (r + 1) 



(/.,. copies) Pr — — • (r + 1) 



(r + 2) 



lr + 2 
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The result: 

After the step t — 1, we have obtained the set Vt-i consisting of representa- 
tions of the form Ljj, j < t, and the quiver representation Mt-i (see fl25|) ). 
which may be considered as a set of representations of the form Lj^. Define 
the representation 

L(A)= U,. (27) 

UjeVt^iUMt-i 

The following proposition will be proved in Section [71 

Proposition 4.1. The representation L,{A) is the canonical form {see The- 
orem of a matrix representation A of the quiver fll2p . 

5 Cycles of linear mappings 

In this section, we give an algorithm for constructing a regularizing decom- 
position that involves only unitary transformations. In the same way, 
one may construct a regularizing decomposition over an arbitrary field using 
elementary transformations. 

By analogy with Section [H we say that a matrix representation A = 
(^1, . . . , At) of a cycle C (see f lT^ ) is regular if 

dimi(A) = ■ • ■ = dimt(A) 

and all the matrices Ai, . . . ,At are nonsingular; otherwise the representation 
is singular. A decomposition 

A~D©---©G©P (28) 

is a regularizing decomposition of A if G are matrix representations 

of the form Gij (see Lemma 13.11) and P is a regular representation. By 
Theorem 13.21 the regularizing decomposition (128!) is determined uniquely up 
to isomorphism of summands. 

The algorithm works like a jack-plane in a woodworker's hands. Starting 
from the vertex 1, we cut a shave: 



16 



We make a full circle by the jack-plane and continue the process until the 
shave breaks away. Then we transpose all matrices of the remaining rep- 
resentation and repeat this process. The obtained representation P of C is 
regular, and the shaves split into a direct sum of matrix representations of 
the form Gjj. 

Note that this proves Theorem 13.21 since P is isomorphic to a matrix 
representation (/„,...,/„, J), where J is a nonsingular Jordan (or Frobenius) 
canonical matrix with respect to similarity; see the end of Section [TJ Hence 
A is isomorphic to a direct sum of representations of the form (i)-(iii) from 
Theorem 13. 2[ The uniqueness of this decomposition follows from Theorem 

The algorithm for cycles: 

This algorithm for every matrix representation 




(29) 



At 

of a cycle C (see ( TT^ ) constructs a decomposition 

A ~ P(A') © A, (30) 

where A' is formed by the matrices of a chain of linear mappings, P sends 
A' to a representation of C that is isomorphic to a direct sum of representa- 
tions of the form Gj^ (see f[T7l) and compare with Example 13.11) . and A is a 
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representation of C that satisfies the following condition for each arrow: 



If the arrow is oriented clockwise, then the matrix 
assigned to it has linearly independent rows. 



(31) 



In steps 1,2,... of the algorithm we will construct quiver representations 

Steps 1, 2, — 1: In step 1 of the algorithm, we check the condition 0311) 
for the representation A and the arrow ai. If this condition holds, we put 
A^^) = A. If this condition holds for a2 too, we put A^^^ = A, and so on. 

If after t steps we found that this condition holds for all arrows of C, then 
we put 

l = t + l, A' = 0, A = A (32) 
and stop the algorithm. Otherwise, we set 



minjz G {1, . . . , t} in A does not satisfy (I3T1) 



(33) 



and continue the algorithm as follows: 



Step I: By unitary transformations at the vertex [/ + 1], we reduce the matrix 
Ai of A to a matrix 







A 



(0 



}d 
}d 



(0 

['+1] 



, rows 



rows 



where the rows of Af^ are linearly independent. This changes we 
denote the obtained matrix by A^''^ and construct the representation 



A«: 



i(0 



(1 + 1)' 



[^ + 1] 




[/ + 2]. 



(the other matrices are the same as in (!29|) ). Its dimensions at the vertices 
{I + ly and [/ + 1] are c^l^+iy and c^l'^^p and the arrow (/ + 1)' — [/ + 2] has 
the orientation of [/ + 1] — + 2]. The matrix ^4^'^ is partitioned into the 

strips ^(j+i)/ and ^[^^ip which are assigned to the arrows (/ + 1)' [I + 2] 

and [/ + 1] — [/ + 2], see (HI). 
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Step r (r > I): Assume we have constructed in step r — 1 a representation 




(r-l) 




.(•■-l) .(---1) 

r I 1 1 ['■+11 r I oi [''+21 
r + 1 r + 2 - 



i('-i) 



where each arrow a^/ : i' {i + 1)' has the orientation of : [i] [i + 1] 

in C, and a^' '■ f' — V + 1] orientation of a\r] '■ [r] — [r + 1]. 

We will reduce A*^*"^^) by unitary transformations at the vertex [r + 1]: 

(i) If a[r] is oriented clockwise, then v4^^'~^^ consists of two vertical strips 
with dim,./ A^^'"-*^) and dimj^j A'^'""^) columns (see f lTS]) - fl2Ul) ): we reduce 
it by unitary row-transformations as follows: 



ir-l)\ /I ('■-1)1 



4(^0 







(34) 



where Af] has linearly independent rows. 



(ii) If a^r] is oriented counterclockwise, then A'^'^~^^ consists of two hori- 
zontal strips with dim^/ A^^"-*^^ and dimj^j A^*""^) rows; we reduce it by 
unitary column-transformations as follows: 



^(r-i) 



\ 1 









L [r] J 


t— > 


* 


At} 



(35) 



where A),, has linearly independent columns. 

3 unitary transformations at the vertex [r + 1] ( 
too; we denote the obtained matrix by A^'^'^ and construct the representation 



These unitary transformations at the vertex [r + 1] change the matrix 
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(/ + 1)' ^ {I + 2)' ^ 



(kt + iy- 

1: 



r'^—{r + iy 



i(r-i) 

[r-l] 



[r]^^[r+l] [r + 2]- 



lr+2] 



-7TT-. {2ty 



(36) 



where A^^^ is partitioned into two strips: 



Air) 



A 



(r) 

(r+1)' 



A 



(r) 
-+1] 



^(r+1)' 



r+1] 



if a[r+i] is oriented clockwise, 



if Q![r+i] is oriented counterclockwise, 



(37) 



and these strips are assigned to the arrows 

(r + 1)' — [r + 2], [r + 1] — [r + 2]. 

The result: 

We make at least t steps and stop at the first representation A^"^ with 



n ^ t and A 



(n) 

(n+iy 



0. 



(38) 



The matrix is assigned to the arrow (n + 1)' — [n + 2]. Deleting 

this arrow, we break A*^") into two representations: 
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A': 



{kt + iy 



{n + iy 



{2ty (39) 



and 



An) 



A: 1 



n 



[?7. + IJ [n + 2\ 



t (40) 



The representation A' is a representation of the quiver 
(^ + 1) 



^{n+iy, (41) 



whose arrows i' (^ + 1)' have the orientation of the arrows : [i] [^+1] 

in C. By analogy with Example 13.11 we construct the mapping P that sends 
a representation B of the quiver (l4Tll to a representation D of the cycle C: 



B 



P 



D 



{i+iy a+2)' ■ ■ ■ — — t' 




^{t+iy 



='(t+2)' 



(fci+1)' -1— (fct+2)'-?iM+2)' ^ („+i)' 



B. 



(fct+i)' 




(42) 



This mapping is known in representation theory as a push-down functor (see 



21 



[H Sect. 14.3]) and is determined as follows: 



I = 



1,2,... 



t 



(43) 



(i.e., Di is the direct sum of all Bji disposed over it), where 



Bi' = Opo with p = dim(;+i)/ B 



(44) 



(recall that the arrow ai is oriented clockwise, see step / of the algorithm), 
and 



(The definition of P : B i-^ D becomes clearer if the representations B and 
D are given by vector spaces and linear mappings: each vector space of D 
is the direct sum of the vector spaces of B disposed over it, and each linear 
mapping of D is determined by the linear mappings of B disposed over it.) 

The following proposition will be proved in Section [HI 

Proposition 5.1. Let the algorithm for circles transform a matrix represen- 
tation A of a cycle C to A' and A. Then 

(a) The condition flHTl) holds for A and all arrows. 

(b) // an arrow ai is oriented counterclockwise and the columns of Ai are 
linearly independent, then the columns of Ai are linearly independent 
too. 

(c) A ~ P(A') © A. 

6 Main theorem 

Theorem 6.1. A regularizing decomposition fl28p of a matrix representation 
A of a cycle C can be constructed in 3 steps using only unitary transforma- 
tions: 

1. Applying the algorithm for cycles to A, we get A c:^ P(A') © A. 
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2. Applying the algorithm for cycles to the matrix representation B := A'^ 
of the cycle {see we get ^ P(B') © B. 

3. Applying the algorithm for chains to A' and M''^ , we get 

A' ~ © • • ■ © L,^,^ , B'^ ~ © • ■ ■ © ■ (45) 

Then ^ 

A ~ Gi^j^ © ■ ■ • © Gi^j^ © B^ (46) 

{see (ITTl) ) and the representation M'^ is regular. 

Proof. By f ITT]) and Proposition 15. 1( c) . 

A ~ P(A') © (P(B'))^ © B^ = P(A') © P(B'^) © B^. 
Substituting we obtain 

A ~ r{U,n) © ■ ■ ■ © P(HiJ ® 

This proves (gG]) sinc^P(Lij) = Gi^-. 

Let us prove that B-^ is regular. By Proposition 15. 1( a) . every matrix of A 
assigned to an arrow oriented clockwise has linearly independent rows. The 
matrix representation B = A^ is constructed by transposing all matrices, 
and it is a representation of the cycle obtained from C by changing the di- 
rection of each arrow. Hence every matrix of B assigned to an arrow oriented 
counterclockwise has linearly independent columns; by Proposition 15.1( b) 
the same holds for the matrices of B. Moreover, by Proposition 15.1( a) every 
matrix of B assigned to an arrow oriented clockwise has linearly independent 
rows. Hence, _ _ _ 

dim[j_|_i] B = rank Bi ^ dim^ B 

for all vertices i = 1, . . . ,t. We have 

dimi B ^ dim2 B ^ ■ • ■ ^ dimt B ^ dimi B. 

Therefore, each matrix Bi is square and its rows or columns are linearly 
independent. So Bi is nonsingular and the representation B is regular. Then 
B"^ is regular too. □ 



23 



7 Proof of Proposition 14.1 



In each step r G {1, ... ,t — 1} of the algorithm for chains (Section H]) we 
constructed the matrix Br of the form (!24|) or fl26l) . Denote by the matrix 
obtained from Br by replacement of all blocks Hi by //. and all blocks * by 
0. Let us prove that the representation 



(r + 1) ^ (r + 2) 



1 



Si 



t 



of the quiver (fT2l) is isomorphic to the initial representation A: 

A~D^, r = l,...,t-l. (47) 
In step 1 we reduced A to 



by unitary transformations at vertices 1 and 2 (see fl22|) ). Using transforma- 
tions at vertex 1, we reduce B\ to 



h 




(48) 



and so A is isomorphic to 
Di : 



1 — - 2 — ^ 3 



A., 



t. 



We may produce at vertex 2 of ©i every transformation given by a nonsin- 
gular block-triangular matrix 



^2 



S\\ 5*12 

^22 



where is k-hj-k if ai : 1 — > 2, and 6*22 is k-hj-k if ai : 1 < — 2. This 
transformation spoils the block Ik of Di but we recover it by transformations 
at vertex 1. 

Reasoning by induction on r, we assume that A is isomorphic to 



Dr-1 K 

r 



(r + ll 



A. 



r+l 



At-i 



t 
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(where Di, . . . , Dj—i are obtained from Bi, . . . , Bj.-i and A'^ is taken from 
(125]) ) and that transformations at vertices 1, . . . , r — 1 may recover the ma- 
trices Di, . . . , Dr-i of Dr_i after each transformation at vertex r given by a 
nonsingular block-triangular matrix 

'S'll 5*12 ■ ■ ■ Sir 
S22 ■ ■ ■ S2r 

Sfp 

in which the sizes of diagonal blocks coincide with the sizes of horizontal 
strips of Br-i if ctr-i '■ ~ 1) — ^ ^, OT with the sizes of vertical strips of 
Br-i if ar-i : (r - 1) < — r, see (124]) and (|26|). 

In step r of the algorithm, we reduced A'^ to S,. of the form or fl2^ by 
unitary transformations at the vertices r and r + 1; moreover, we used only 
those transformation at vertex r that were given by unitary block- diagonal 
matrices partitioned as (H9!) . By the same transformations at the vertices 
r and r + 1 of Dr-i, we reduce its matrix A'^ to Br. Then we reduce Br 
to Dr by a transformation at vertex r given by a matrix of the form fl49p . 
and restore Di, . . . , Dr-i by transformations at vertices 1, . . . , r — 1. The 
obtained representation is D^, and so 

A ~ ^r-l ^ Or- 

Moreover, we may produce at the vertex r + 1 of all transformations given 
by block-triangular matrices, restoring the matrix Dr by transformations at 
vertex r given by matrices of the form fH9l) . and then restoring Di, . . . , Dr-i 
by transformations at the vertices 1, . . . , r — 1. This proves the isomorphism 

We now transform the representation (1 ^ f < ^) of the quiver f|T2|) to 
a representation of a new quiver as follows. We first replace each vertex 
i G {1, 2, . . . , r + 1} of by the vertices ii, . . . , id^, where 

di = dimj = dim^ A. 

A' 

Then we replace the arrow (r + 1) — ^ (7- _|_ 2) by the arrows 



(49) 
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(r + l)i 



(r + 2) 



of the same direction. 

At last, we replace each arrow i — — (i + I) with i ^ r hj arrows that 
are in one-to-one correspondence with the units of the matrix Di: every unit 
at the place [p, q) in Di determines the arrow 

iq^{i + l)p if di-i — ^(2 + 1), 



or the arrow 



ip^{i + l)q if ai:i< — (^ + 1)- 



(These arrows represent the action on the basic vectors of the linear operator 
Oi © ■ ■ ■ © Cid, C(z + l)i © • ■ ■ © C(2 + 

directed as : i — + and given by the matrix Di.) Since in each row 
and in each column of Di at most one entry is 1 and the others are 0, two 

arrows ip (z + l)g and y (z + 1),' have no common vertices {Vi sends 

each basic vector to a basic vector or to and cannot send two basic vectors 
to the same basic vector). Denote the obtained representation by Q,.. 

The quiver representation Qt_i is a union of nonintersecting chains; each 
of them determines a representation of the form Ljj . Hence, is a direct 
sum of these representations. By fH7|) . A O^-i, so Dj^i is the canonical 
form of A, and we need to prove L(A) = Dt_i (see fl7r|) ). 

It suffices to show that 



Qt-i = Vt-i U Mt-i 

(see the set of indices in (1271)). The equality Qi = Pi U Mi holds since the 
matrix Di is obtained from Bi by replacement of H with Ik (see (l22l) and 
( HHj) ). Reasoning by induction, we assume that Qr-i = Pr-i U M^-i. Then 
= U Mr by the construction of Vr and in step r of the algorithm 
for chains and since D^ is obtained from by replacement of all blocks Hi 
by and all blocks * by 0. This proves Proposition 14. 1[ 
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Example 7.1. Suppose we apply the algorithm to a matrix representation 



of dimension (4,5,4,5) and obtain 

^ ^ O31 Hi 
O21 O23 



where Hi, . . . , are nonsingular 3 x 3, 2 x 2, 1 x 1, 2 x 2, and 1x1 matrices. 
Then 





' 021 


H2 


* 


* 




' O22 




* 




On 


O12 


On 


^3 


, B,= 


O12 


O12 


* 




On 


O12 


On 


On 




O12 


O12 





: 



10 

10 

1 







10 

10 

1 





10 

1 



1 



and 



ll 


^^2i 


^^3i 


4i 




^^22-^ 


^.32 


42 


Q3 : I3 ^^^'^^ 


^^23-^ 


33 


43 




24 


34^ 


44 
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^^45 


We have the canonical form of A: 






A =i Ln © IL12 i 


D Li4 © Li4 i 


D L22 © L23 © L34 i 


D L44 © L44 



(50) 



Note that the block-triangular form of 5"^ (see 
disposition of the chains 



follows from the 



h h 
Pi {Pi + 1) 



,r + l, 
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in the quiver representation M^-i (see f l23l) ): they represent the hnear map- 
pings and we may add these chains from the top down by changing bases in 
vector spaces; this is clear for the quiver representation flSO]) . 



8 Proof of Proposition 15.1 

The representation A*^*'-' (see f l36|) ) is a representation of the quiver, which 
we will denote by Q^**-*. For every representation 




of this quiver, we define the representation 

F(B): 1 — 2 t 



Dt 

of the cycle C by "gluing down of the shave" (see the beginning of Section 
ED: 

^Bi ifi^[r+l]. 



B ifi=[r+l]. 



where Bi^ is defined by (jH]) (compare with 0431) ) . The mapping F is analogous 
to the "push-down functor" fH21) . Moreover, for the representation A'^"\ 
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obtained in the last step of the algorithm for cycles, we have 

F(A(")) =P(A') ©A, 

where A' and A are the representations fl5U]) and fHU]) . 
By (l20l) . the matrix B in (!5T!) has the form 



(52) 



B(^r+iy I -B[r+i] ] if «[r+i] is oriented clockwise, 



5 



(r+l)' 



5^. 



r+l] 



if <y\r+i] is oriented counterclockwise. 



By triangular transformations with a representation B of the form (15T]) . we 
mean the following transformations: 

(i) additions of linear combinations of columns of -B[r+i] to columns of 
B(^r+iy if tt[r+i] is oriented clockwise, 

(ii) additions of linear combinations of rows of -B(r+i)' to rows of -B[r+i] if 
a[r+i] is oriented counterclockwise. 

We say that B is a triangular representation if 

F(B) ~ F(B^) 

for every representation B^ obtained from B by triangular transformations. 

Lemma 8.1. Suppose D is obtained from a triangular representation B of 
Q^^'^ by transformations at the vertex [r + 2]. Then D is triangular too. 



Proof. Let 

5 = (J,...,/, 5[,+2], /,...,/) :B^D 

(see dl])). We must prove that F(D) ~ F(D^) for every obtained from 
D by triangular transformations. Denote by B^ the matrix representation 
obtained from B by the same triangular transformations. By ([8]) and the 
definition of triangular transformations, there is a block matrix 



R 



I 
* I 
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such that 



[^(^■+1)' I ^[r+i]] = S[r+2][B(r+iy \ B[r+i]]R if ^[r+i] IS orieiited clockwise, 

SZ^2] if ci[r+i] is oriented counterclockwise, 




B[r+2]S^.^.^2] if is oricutcd clockwise, 

S[r+2]Bir+2] if is orieutcd counterclockwise. 



These equahties imply 

5=(/,...,/,5[,+2],/,...,/):B^^D^ 

and 

F(D) ~ F(B) ~ F(B^) ~ F(D^). 

□ 

Lemma 8.2. Each representation A^'''^ {obtained in step r of the algorithm 
for cycles) is triangular and F(A*^^''') A. 

Proof. The lemma is obvious if I = t + 1 (see (|32l) ). Suppose / ^ t. The 
statements hold for A^^\ . . . , A'^'\ Reasoning by induction, we assume that 
they hold for A*^^'"^^ with r — 1 ^ I and prove them for A*^*"^. 

First we apply the unitary transformations at the vertex [r + 1] from step 
r of the algorithm for cycles to the representation A*^^~^) of the quiver Q^'""^): 
we reduce the matrix A^^'^^^ to a block-triangular form by transformations 
or ([35]) (depending on the orientation of a[r]), and the matrix ^[^+1]^ to 
A^^\ Denote the obtained representation by A^'^^'^^^\ 

Then we make zero the block * of fl34|) or fl35|) by triangular transfor- 
mations and obtain the following representation A^^'^^^^^ of the quiver Q^^~^^: 



R 



(r+l)' 



B[H 
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By the induction hypothesis, A ~ F(A*^'" and A*^^' is triangular. By 
Lemma I8.H A*^*""^/^) is triangular too, and so 

F(A("-2/3)) ^ F(A("-i/=^)). 

We have 

A ~ F(A(''"^)) ~ F(A(''~2/3)^) _ f(A(^''"^/^)) = F(A('^)). 

Let A*^^-*^ be obtained from A*^**) by triangular transformations. These 
transformations reduce A^^^ (see ( l36l) ) to a new matrix and do not 

change the other matrices of A^^^^. Since 

Air) _ A'-yS) 
^ - ^[r+l] ' 

these transformations with can be realized by transformations at the 

vertex [r + 1] of A*^''"^/'^-'; denote the obtained representation by A^'^~^^^'>^, it is 
triangular by Lemma [8. 1[ These transformations may spoil the subdiagonal 
block of 

4 ('^-1/3) _ Ai"-) ff. 4W 
^ — ^r' ^ ^[r] ' 

but it is recovered by triangular transformations and so 
Since 

f(a(")) = f(a('-^/^)) ~ f(a("-i/3)^) ~ F(A(")^), 

the representation F(A*^^''') is triangular. □ 
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Lemma 8.3. Let A^''^ be the representation obtained from a representation 
A in step k of the algorithm for cycles, and let k ^ I {hence I ^ t by 
and fl33|) ). Denote 

.(,)^jAf) tft^[k + l], 
if i= [k + 1], 



A. 



where i = 1, . . . ,t. Then 

(i) The rows of Af'^ are linearly independent if cXi is oriented clockwise and 
i ^ k. 

(ii) The columns of Af^"* are linearly independent if ai is oriented counter- 
clockwise and the columns of Ai are linearly independent. 

Proof. We will prove the lemma by induction on k. Clearly, the statements 
(i) and (ii) hold for k = I. Assume they hold for k = r — 1 ^ I and prove 
them for k = r. We need to check (i) and (ii) only for i = [r] and i = [r + 1] 
since in step r of the algorithm we change A^^j and A^^'^^y . 

By the matrix A^^j = A^^^ has linearly independent rows if a[r] is 
oriented clockwise. By fl5S]) . this matrix has linearly independent columns 
if both a[r] is oriented counterclockwise and ^[r]"^'' = A^^'^^^ has linearly 
independent columns. Hence, (i) and (ii) hold for i = [r]. 

The statements (i) and (ii) hold for i = [r + 1] by the induction hypothesis 
and since ^[^+1] = ^^^^ is obtained from ^[^+1]'' elementary transformations 
with its columns or rows. □ 

Proof of Proposition 15. it The statement (c) of Proposition 15.11 follows from 
fl52|) and Lemma W2[ so we will prove (a) and (b). 

If / = t + 1 (see ([32])),^hen A = A satisfies (a) and (b). 

Suppose / ^ t. Then A is the restriction of the representation A*^") (ob- 
tained in the last step of the algorithm) to the cycle C and so Ai = A[^'' 
(i = l,2,...,t). 

Since 

^{n) ^ ^(n) ^ ^ 



if i ^ [n+ 1], 



T(n) _ 4(„,) 



An) 



^[n+1] 
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if «[„+!] is oriented clockwise (see ( 1571) and (1551) ). and 



7in) _ Un) 








■ " 


. ["+11 . 







if a[n+i] is oriented counterclockwise, the statements (i) and (ii) follow from 



Lemma 18.31 in which k = n ^ t. 



□ 



The author wishes to express his gratitude to Professor Roger Horn for 
the hospitality and stimulating discussions. 



References 

[1] H. Bass, Algebraic K-theory, Benjamin-Addison Wesley, 1968. 

[2] G. Belitskii, Normal forms in matrix spaces. Integral Equations and Op- 
erator Theory, 38 (no. 3) (2000) 251-283. 

[3] G. R. Belitskii and V. V. Sergeichuk, Complexity of matrix problems. 
Linear Algebra AppL, 361 (2003) 203-222. 

[4] N. M. Dobrovol'skaya and V. A. Ponomarev, A pair of counter-operators 
(in Russian), Uspehi Mat. Nauk20 (no. 6) (1965) 80-86; MR 36#2631. 

[5] P. Donovan and M. R. Freislich, The Representation Theory of Finite 
Graphs and Associated Algebras, Carleton Lecture Notes 5, Ottawa, 
1973. 

[6] P. Gabriel, Unzerlegbare Darstellungen I, Manuscripta Math. 6 (1972) 
71-103. 

[7] P. Gabriel and A. V. Roiter, Representations of Finite- Dimensional Al- 
gebras, Encyclopaedia of Math. Sci., Vol 73 (Algebra VIII), Springer, 
1992. 

[8] F. R. Gantmacher, The Theory of Matrices, Vol. 1, Chelsea, New York, 
1971. 

[9] J. Gelonch, The contragredient equivalence for several matrices: a set 
of invariants. Linear Algebra AppL, 291 (1999) 37-49. 



33 



[10] O. Holtz, Applications of the duality method to generalizations of the 
Jordan canonical form, Linear Algebra AppL, 310 (2000) 11-17. 

[11] R. A. Horn and D. I. Merino, Contragredient equivalence: a canonical 
form and some apphcations. Linear Algebra Appl. 214 (1995) 43-92. 

[12] L. Kronecker, Algebraische Reduktion der Scharen bilinearer Formen, 
Sitzungsber. Akademie Berlin (1890) 763-776. 

[13] L. A. Nazarova, Integer representations of a four-group (in Russian), 
Dokl. Akad. Nauk SSSR 140 (no. 5) (1961) 1101-1014; MR 24#A770. 

[14] L. A. Nazarova, Representations of a tetrad (in Russian), Izv. Akad. 
Nauk SSSR Ser. Mat. 31 (no. 6) (1967) 1361-1378; MR 36#6400. 

[15] L. A. Nazarova, Representations of quivers of infinite type. Math. USSR 
Izv. 7 (1973) 749-792. 

[16] L. A. Nazarova and A. V. Roiter, Finitely generated modules over a 
dyad of two local Dedekind rings, and finite groups which possess an 
abelian normal divisor of index p (in Russian), Izv. Akad. Nauk SSSR 
Ser. Mat. 33 (1969) 65-89; MR 41#5479. 

[17] P. Rubio and J. Gelonch, Doubly multiplicable matrices. Rend. Inst. 
Mat. Univ. Trieste 24 (no. 2) (1992) 103-126. 

[18] V. V. Scrgeichuk, Canonical matrices for linear matrix problems. Linear 
Algebra Appl. 317 (2000) 53-102. 

[19] P. Van Dooren, The computation of Kronecker's canonical form of a 
singular pencil. Linear Algebra Appl. 27 (1979) 103-140. 



34 



